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From pressure and surface dilation measurements, we show that a sohd-hquid-type transition 
occurs at low excitation frequencies in vertically vibrated granular layers. This transition precedes 
subharmonic bifurcations from flat surface to standing wave patterns, indicating that these waves are 
in fact associated with the fluid like behavior of the layer. In the limit of high excitation frequencies, 
we show that a new kind of subharmonic waves can be distinguished. These waves do not involve 
any lateral transfer of grains within the layer and correspond to excitations for which the layer 
slightly bends alternately in time and space. These bending waves have very low amplitude and we 
observe them in a vibrated two-dimensional layer of photoelastic particles. 
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I. INTRODUCTION 

In a driven granular system, such as a vibrated layer 
composed of macroscopic particles, energy is dissipated 
by inelastic grains collisions and by friction. Also, for 
realistic situations, the typical energy of a grain is many 
orders of magnitude greater than ksT, so the temper- 
ature does not play an important role in the dynamics 
of these materials. However, it is commonly observed 
that granular matter behaves like solids, liquids or even 
gases depending on the energy injection and energy dis- 
sipation rates 0|. More precisely, the competition be- 
tween these quantities determines the residual velocity 
fluctuations which in turn play the role of thermal fluc- 
tuations in these materials. From this view point, sur- 
face waves excited by vertical vibrations provide one of 
the most striking examples in which granular materials 
act like a fluid. For a granular layer, it was well estab- 
lished that parametric waves can be observed whenever 
the dimensionless acceleration F = A(27r/)^/f/ exceeds a 
critical value |^ (here g is the acceleration of gravity, A 
is the amplitude of the vibrating surface, and / is the 
frequency of the driving force). The primary instabil- 
ity resulted from a flat layer to a pattern of squares or 
stripes, depending on both / and the particle diameter, 
d. As an illustration we present two typical snapshots of 
these kind of waves in Fig. |l] a and b. It was found that 
the crossover frequency fd at which the square to stripe 
transition occurs is proportional to dr^^"^. In addition, 
this scaling was qualitatively understood in terms of the 
ratio of kinetic energy injected into the layer to the po- 
tential energy required to raise a particle a fraction of 
its diameter. At low /, the horizontal mobility should be 
high since the layer dilation is large , whereas at large 
/, the mobility should be low since layer dilation is small 
We notice, however, that these waves are the corre- 
sponding hydrodynamic surface modes of the layer since 
they require a transfer of mass to be sustained |^,^,^. 

In this paper, we report on measurements of the pres- 



sure due to the layer-container collision and the surface 
and bulk dilation of the layer. We show that there exists 
a critical F for which the flat layer undergoes a phase 
transition. This critical F is smaller than the critical one 
at which waves appear. At low / this transition is a solid- 
liquid type, at intermediate / only a heating up of the 
surface layer is observed, whereas at high / a compaction 
transition is detected. Thus, our results show that gran- 
ular surface waves are naturally linked to the fluid like 
behavior since they are observed only when the energy in- 
put per particle is enough to induce a minimal dilation of 
the flat layer. At high / and for the same critical value of 
F at which hydrodynamic waves appear, bending waves 
are detected instead. In this regime, our measurements 
show that the mobility of grains is almost completely sup- 
pressed in both the bulk and the free surface of the layer. 
In this case, a set of experiments conducted with large 
photoelastic particles allow us to observe these waves as 
an alternation of bright and dark zones that oscillate in 
time at half of the forcing frequency. In addition, we show 
for the low frequency regime that when F is increased to 
a value close to 4.6, an inverse transition from a fluid to a 
compact layer is observed. This transition is responsible 
for the flat with kinks state reported in previous works 
Q (See Fig. |^ c and d). A brief summary of our results 
has been given in ref. H. 

This article is organized as follows: section II is de- 
voted to the description of our experimental setup. In 
section III, which is the main part of this article, we 
present our pressure measurements and we establish the 
basis to obtain density profiles and the bulk dilation from 
the experimental data. In section IV we present surface 
dilation measurements and we discuss how they link to 
the bulk dilation measurements. The energy dissipation 
rate in the vibrated layer is also estimated as a function 
of the excitation frequency. Section V is concerned with 
the study of the very low amplitude regime of surface 
waves and its connection with the compact state of the 
layer. Finally a brief discussion is given. 
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FIG. 1. Snapshots of the layer state at differents conditions obtained in a large cell, a) Squares, / 20 Hz, F — 3.0. b) 
Stripes, / ~ 60 Hz, F = 3.0. c) Flat with kinks at / ~ 25 Hz, F = 4.7. d) Flat with kinks at / ~ 60 Hz, F = 4.7. 



II. EXPERIMENTAL SETUP 

In this article we present two sets of experiments. In 
the first set, we measure simultaneously the pressure re- 
sulting from the collisions of the layer with the vibrating 
plate and the surface particle dilation. From the pressure 
measurements, we obtain information on the bulk dila- 
tion and density (see section III). The surface dilation is 
determined through the normalized reflected intensity of 
the surface layer (See discussion below). 

In fig. ^ we present a schematic drawing of the ex- 
perimental setup. In this experiment, a thin layer of 
0.106 — 0.125 mm diameter bronze particles, 15 particles 
deep, is placed at the bottom of a 40 mm diameter and 25 
mm height cylindrical container. The container's wall is 
Lucite while the base is aluminum to reduce electrostatic 
effects. The container is mounted on a high frequency 
response pressure sensor (PCB. Model 208A11) which is 
driven by an electromechanical vibration exciter. The re- 
sulting acceleration is measured to a resolution of O.Olg. 
A second Lucite cylinder is used as a lid for the whole 
system allowing evacuation of the container to less than 
0.1 Torr; at this value volumetric effects of the gas are 
negligible [p|. The surface of the layer is illuminated at 
low angle (20^ respect to the horizontal) by an array of 
18 LEDs organized in a 10 cm diameter ring. The re- 
flected light from the surface layer is focused by a lens of 
28 mm focal length on a flat photodiodc of 25 mm^ area. 
The whole system is automatically run by a Power PC 
computer equipped with A/D and GPIB boards. 

With this set-up, the measured light is proportional to 
both the incident light and the reflectivity coefhcient of 
the bronze particles i? (i? ~ 0.6). We notice that only a 
small fraction s, which is about 5 percent of the surface of 
a single particle, reflects light in the direction of the solid 
angle of the camera (in the experiment, lens aperture 
angle is about 12°). The intensity /, measured by the 
photodiodc, can be then taken proportional to the sur- 
face density or more precisely to the number of particles 
within the first layer. Furthermore, / relates to surface 
dilation Sg as I/Iq « (P /{d + Sg)"^, where Iq is a refer- 
ence intensity for which we take (5s = [Q. We neglect 
multiple reflections, since their dominant contribution is 
proportional to both and s^, where represents the 
probability of having a secondary reflection within the 



solid angle of the camera. Thus, by taking the incident 
light in a small enough angle, we insure that, to a first 
approximation, only the first layer of particles contribute 
to /. 

The second set of experiments consists of vibrating a 
two dimensional layer of photoelastic particles and tak- 
ing images of the layer motion with a high speed CCD 
camera. The goal is to observe the very low amplitude 
waves found in the high frequency regime. A cell made 
of two glass plates 400 mm wide by 100 mm high was 
mounted on the moving platform of our vibrator system. 
The gap between the plates is controlled by spacers of 
varying thickness to a resolution of 0.05 mm. Images 
with resolution of 256 x 256 pixels are captured at rates 
of 1200 frames per second by a Hisis 2002 CCD camera. 
Images are obtained by transmission using parallel light 
with incidence perpendicular to the cell. Transmitted 
light is filtered through a polarizer whose axis is perpen- 
dicular to direction of polarization of the incident light. 
More details are given in section V. 




Accelerometer 
M Mask 

LED's 
Incident Light 



- Force Sensor 



Accelerometer 



Force Sensor 



' z(l) = A sin(wt) 



FIG. 2. A schematic drawing of the apparatus showing the 
cylindrical cell, the location of the pressure sensor, and the 
setup for reflectivity measurements. 
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III. PRESSURE MEASUREMENTS AND 
DENSITY PROFILES 

A. Experimental procedures and collision model 

In Fig. I^a, we present a typical pressure signal P{t) as 
a function of time for / = 40 Hz and F — 2.2. This signal 
is composed of a sinusoidal component, corresponding to 
the force required to accelerate the cell, and of a peak 
sequence, due to the layer-plate collisions. We are inter- 
ested in the shape of the peak during the collision so it is 
necessary to subtract the sinusoidal component. By this 
procedure we obtain a number of pressure peaks, typi- 
cally about 15. As an illustration, we show in Fig. yb the 
average peak obtained from the peaks presented in Fig. 

To characterize the peaks of a sequence like the one 
presented we take the maximum value P and the width 
Tc, which is a measure of the collision time. More pre- 
cisely these quantities are obtained as the averages of the 
sequence {pW,ri'=^}, where fc is an integer that indicates 
the collision number, i.e. P = {P^''^) and = {^ri''^^ 

where ( ) denotes the sequence average. For the data in 
the Fig. I, P = 2.45 ± 0.03 kPa and = 0.52 ± 0.01 ms 
(In this case Tc is taken as the width at a pressure equal 
to a quarter of P). 
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FIG. 3. (a) Time evolution of the pressure signal P{t) for 
/ = 40 Hz and F = 2.2. In (b) we present the averaged colli- 
sion peak obtained from P{t), and we show the definition of 
the maximum pressure P. The collision time Tc is the width 
of the peak at certain height. 

Now, due to the impulsive nature of the layer plate 
collision, we write the following scaling relation 



Ap Tc 



(1) 



where AI is the total mass of the layer, Ap is the active 
surface of the plate and Vc is the layer-plate relative ve- 
locity at the collision. If the collision takes place at time 
t = i and is completely inelastic then Vc = VpiJ!) — Vbit) , 
where Vp and Vf, are the plate and layer velocity respec- 
tively. Thus, in practice, we approximate the collision 
velocity by the one predicted by the completely inelastic 
ball model. We note that this model assumes that the 



layer losses immediately all its energy and takes the ve- 
locity of the moving plate. This approximation is valid 
since previous experiments have shown, by measuring the 
flight time of the layer, that the motion of the center of 
mass of the layer follows the motion of an inelastic ball 
g|. 

We remark that in our experiments, even close to 
F ~ 1, the vertical average bulk dilation of the layer, 
5h, is never strictly zero. This is concluded because the 
collision time is not determined by the Hertz theory. In- 
deed, experiments performed by E. Falcon et al ^ have 
shown that, in the case of a column of N particles in con- 
tact {5b = 0) colliding vertically with a fixed plate, the 
collision time is Tc — {N — l)Tq + Ti, where N is the num- 
ber of particles, Tq is the time duration of the momentum 
transfer from one particle to another, and ti is the colli- 
sion time between two particles predicted by Hertz theory 
Ip^ . Now, for the particles used in our experiments the 
orders of magnitude of these times are « ri « 10~^ 
s [O. Then, the predicted collision time for our three 



dimensional layer should be of order Tc 



10- 



This 



order of magnitude is never observed in our experiments, 
where for F close to 1 and for a wide range of / the min- 
imum collision time is of order Tc « 5 x 10"'' s. The 
difference between what is expected from the Hertz pre- 
diction for (5f, = and the experimental values for Tc 
can be explained if the layer is slightly dilated at the 
collision; a local dilation of 1 percent {5b/ d ^ 0.01) dras- 
tically changes the collision regime from Hertz contact 
to ballistic collisions. Although this effect has interest- 
ing consequences on sound propagation, it does not affect 
the center of mass motion of the layer. Thus, for F > 1 
we write the total dilation of the layer as A = VcTc. If 
in addition we assume that dilation is homogeneous, i.e. 
A = N5b, where N is the number of layers, we can write 
equation (Q) as 



P 



M If 
'A'plUb 



(2) 



Thus, our procedure to measure the local bulk dilation 5b 
is: we fit P versus Vc/Tc to obtain the numerical factor 
of equation (|l|) and we then use equation (||) to obtain 
the value of 5b- 

In the following, we will relate the layer density to the 
pressure signal P{t). Fig. |^ presents a schematic view 
of a collision. We consider that the plate collides with 
the granular layer of density p{z) at a velocity Vc at time 
t = 0. Then, the layer is initially fixed in space and the 
plane {x,y,z = 0) of the reference frame {x,y,z) coin- 
cides with the inferior layer just before the collision (See 
Fig. IJa). We also notice that the mass of the plate Mp is 
much larger than the mass M of the layer {M/Mp ~ 0.01) 
so we neglect any velocity change of the plate due to mo- 
mentum transfer. We also neglect the variation of the 
relative velocity due to the action of gravity since exper- 
imentally we observe that gTc « Vc- Thus, the force 
exerted on the plate due to de momentum transfer is 
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FIG. 4. Schematic of the colhsion. Fig. (a) shows the in- 
stant t = at which the plate colides with the layer. The 
velocity of the collision is Vc and the density of the layer 
before the collision is p{z)\ Fig. (b) shows how the layer ac- 
cumulates over the moving plate at an instant t > 0. We also 
present the definition of Zc{t) as the position of the compre- 
sion front relative to the origin of the fixed frame {x,y,z). 
In all the calculations the effect of gravity is neglected, since 
experimentally gTc <^ Vc- 



dm{t) 
dt 



Vr. 



(3) 



Here m(t) is the mass that has collided with the plate 
at time t. In the next we need to link dm{t)/dt to p{z), 
measured in the fixed frame {x^y, z). To do this we de- 
fine the compression front Zc{t) as the height (relative 
to the fixed frame) of the mass that has been deposited 
on the plate at time t (see Fig. ||b). So, the mass m(t) 
accumulated on the plate becomes 



i{t) — Ap I p{z)dz 



(4) 



whose derivative with respect to time is 

dm{t) /+^^^^c(^) 
- App{zc[t))- 



dt 



dt 



(5) 



Now, at time t the plate has moved a distance Vet from its 
initial position, and the height of the layer with respect 
to the plate can be written 



h{t) 



m{t) 
ApPo 



— / p(z)dz 

Po J 



(6) 



where po is the density of the layer in the compact state. 
We then obtain 



Zc(t) = Vct + 



1 

Po 



=(*) 



p{z)dz 



(7) 



Differentiating this equation respect to time, we obtain 
the density evaluated at the location of the compression 
front as a function of its velocity Zc(t) 



Zc{t) 



(8) 



Using this result in equation (|^), the relation (||) for the 
force becomes 



m 

ApPoVc 



(9) 



Integrating this equation between and t gives us an ex- 
pression for the compression front as a function of time 



Zc{t) = V.t + 



1 



ApPoVc 



F{t)dt 



(10) 



In summary, from the pressure signal P{t) ~ F(t)/Ap we 
deduce the time evolution of the compression front Zc{t). 
Also, we use equation (^), which tells us about the time 
evolution of density at the compression front. Since the 
time t enters as a simple parameter we can obtain the 
density as a function of Zc, right before the collision. 

Finally, we notice that both equations (^) and (||) are 
laws for the conservation of mass and momentum through 
the compression front. It is possible to deduce from them 
the velocity of the compression front and the pressure in 
the compressed part of the layer as implicit functions of 



Zc{t) 

m 

A, 



PoVc 



Po - p{Zc{t)) 
Pop{Zc{t))V^ 
Po - P{Zc{t)) 



These expressions are generalizations of those obtained 
by Goldshtein et al |12| for the propagation of a shock 
wave through an homogeneous layer of inelastic particles. 

Before presenting our exp erimental results we will 
show that with equation (y^ we can check the momen- 
tum conservation of the collision. Evaluating it at time 
t = Tc (here Tc is the total collision time) we obtain 



1 



ApPoVc 



F{t)dt 



(11) 
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Since by definition Zc{Tc) = H + VcTc, where H is the 
layer thickness in the compact state and VcTc is the to- 
tal displacement of the plate during the collision, and 
Po = M/HAp, we find 

j F{t)dt = MVc (12) 



B. Experimental results 

We begin this section by presenting results concerning 
the momentum conservation relation (^2|). Experimen- 
tally we have checked it for a wide range of parameters, 
namely 1 < T < w 2.8 and 35 Hz < / < 350 Hz. The 
main point is that, in our experiment, the velocity Vc cal- 
culated from the completely inelastic ball model is a very 
good approximation. As the internal degrees of freedom 
of the layer are excited one expects that this approxima- 
tion becomes less accurate. However, this only occurs in 
the wave regime where the take off velocity is reduced 
and then both the collision velocity and the flight time 
are smaller than predicted |^ . 

There are other experimental effects that should be 
considered. For instance, friction on the cell walls can 
tranfer momentum to the layer. Another posibility, much 
less probable, is the transfer of momentum to the walls 
by the formation of dynamical arcs. Therefore, we ex- 
pect that the correct form of the momentum conserva- 
tion, taking into account all possible sources of errors, to 
be 



Jc 

J F{t)dt = M.ffV, 



(13) 



where Me// is an effective mass, and Vc is the velocity 
of collision calculated from the completely inelastic ball 
model. From this model, we know that Vc = gF{T)/f, 
were F{r) is a nonanalytic function of F, which is found 
numerically [^. From our data we find then that the 

Tc 

scaling / F{t)dt ~ Vc is very well verified as a function 



of both F and /. As expected, we find that Afg// is 
independent of both F and / and is slightly lower than 
M. 

In order to complete our description, it is necessary to 
estimate the granular density in the compact state po- 
This quantity is in principle a dynamical variable, in the 
sense that it depends on how we compact our layer. How- 
ever, we notice that the important parameter in ( p^ is 
ApPo = M/H. Taking in to account the previous dis- 
cussion about the effective mass of the layer and that 
H X 1.7 mm, we obtain ApPo ~ 4.2 kg/m. 



Fig. H presents averaged pressure peaks for several 
values of F at / = 40 Hz. Each curve is obtained by av- 
eraging 15 collisions, in the same way as that for the one 
presented in Fig. |3|3. We notice that our analysis is valid 
for cases for which F < F^, where w 2.8 is the onset 
for parametric waves. Thus, the last pressure curve pre- 
sented for F = 2.98 is shown to display the difference of 
the parametric wave state; at this value of F the collision 
is quite spread out in time and the maximum pressure 
achieved is also much smaller. We will discuss more this 
point in the next section. For the other pressure curves 
presented, the intensity of the collision seems to be an 
increasing function of F, i.e. the maximum pressure P 
increases. We understand this as the simple fact that the 
relative velocity of the collision Vc is increasing with F, as 
it does in this region of F in the completely inelastic ball 
model. What is not possible to explain with this model 
is the observation that the collision time Tc also seems to 
be an increasing function of F. As discussed before, this 
is due to the excitation of the internal degrees of freedom 
of the granular layer, i.e, the layer is dilated. 

In Fig. ^ we present the layer density as a function 
of height for each of the curves introduced in Fig. ||. In 
general, the density p{z) is approximately constant in a 
center region and decreases toward both ends of the layer. 
A penetration length for the dilation of the layer can be 
identified. This length increases with F; for instance, it 
is of the order of d and 3d for F « 1.4 and F « 2.4 re- 
spectively. Another fact that traces back the dilatance of 
the layer is that its thickness increases with F; it is 14(i 
for F « 1.4 and 17d for F « 2.8. In the limit of F - 1 we 
obtain that the density in the central part of the layer is 
close to Po- 

Another interesting piece of information is the bulk 
dilation within the layer, S, as a function of z. This 
quantity is linked directly to the density as p{z) — 
i>rao/(d -\- (5)^, where v is an average coordination num- 
ber and rrio the particle mass ||]. Defining (5 = for 
p{z) — Po we obtain 



d 



P(f) 

Po 



- 1 



(14) 



Then , it is possible to find S{z) for all the data pre- 
sented previously, which is what is shown in Fig. |^. 
From this figure, we can qualitatively examine the de- 
pendence of 6(z) on F, as we did with p{z). For instance, 
as we increase F, the top and the bottom of the layer 
continuously dilate and the extension of the dilated part 
increases. The central part also dilates as F increases; 
for F « 1.4 and F « 2.4 it is of the order of 0.004d and 
0.02d respectively. We also notice that for F « 2.4, the 
value 6 w O.ld is reached for z k. 2d and z k, 15(i, while 
the total height of the layer is approximately lid. 
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FIG. 5. Averaged pressure peaks as functions of time for various values of T and / = 40 Hz. For F = 2.98 the layer is already 
in the wave regime. 




FIG. 6. Normalized density, p{z)/po, as function of height, z, for various values of F and / = 40 Hz. 
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FIG. 7. Normalized dilation, S{z)/d, as function of height, z, for various values of F and / = 40 Hz. 



At this point, it seems appropriate to define a quan- 
tity as the vertical average of S{z); we do this because it 
is necessary to hnk it in some way to Sf,. Therefore, we 
define this average dilation as 



H' 



(15) 



where H' is the total height of the layer. Notice that H' 
depends on the excitation intensity (i.e, it depends on F 
and /). 

Fig. ^ presents {6) versus F for several excitation fre- 
quencies. We observe in the low frequency regime that 
{S) presents a transition at F « 2. In contrast, at high 
/ this transition has the tendency to disappear and we 
observe that (S) is roughly a linear function of F. Com- 
plementing the previous data, in the inset of Fig. ^ we 
present {6) versus / for two values of F, below and above 
F = 2. The average dilation {6) is a decreasing func- 
tion of /. We also present with continuous lines the fits 
{6) ~ we find b = -1.38 ± 0.06 and b = -1.58 ± 0.03 
for F = 1.5 and F = 2.4 respectively. Similar behaviors 
are obtained for Sf,, as a function of both F and / (See 
section IV). 

To conclude, in this section we have introduced the 
peak shape of pressure from which we have obtained den- 
sity profiles. We have discussed these results in a qual- 
itative way and found that dilation is a function of the 
vertical coordinate. These results, which show that dila- 
tion is approximately constant in the bulk but increases 
sharply close to the free surface of the layer, are similar 



to the ones reported in previous works in one and two 
dimensions [|l3|,0. We have also shown that at a low 
frequency of excitation, a critical value of F w 2 exists 
where an abrupt change in the layer dilation takes place. 
This transition will be discussed in detail in the next sec- 
tion. 
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FIG. 8. Average dilation, (S), versus F for various values 
of /; (.) / = 40 Hz, (A) / = 55 Hz, (+) / = 70 Hz, (o) 
/ = 120 Hz and (*) / = 250 Hz. In the inset, log-log plot 
of (5) versus / for F = 1.5 (•) and 2.4 (A). The continuous 
lines show the fits (S) ~ f for each F, with b = -1.38 ± 0.06 
and b — —1.58 ± 0.03 respectively. 
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IV. PRESSURE AND REFLECTIVITY 
MEASUREMENTS 

Both the time-evolution of the pressure and the inten- 
sity are represented in Fig. ^ as a function of F when / is 
in the low frequency regime. For F > 1, the layer-plate 
collision is always visible in the pressure signal as large 
peaks. However, no trace of this collision is observed in 
the reflected light up to F - 2. Thus, for 1 < F < 2 
the layer is compact, indicating that the energy injected 
during the layer-plate collision is completely dissipated 
by the multiple collisions between the grains or by fric- 
tion. In contrast, for F > 2 the time mean value of the 
reflected light (DC component) exhibits a strong decrease 
that shows that the layer undergoes a transition from a 
compact to a dilated state. A modulation in time (AC 
component) which oscillates at the forcing frequency is 
also observed in the reflected light for F > 2. This mod- 
ulation is in phase with the pressure peak. Indeed, im- 
mediately after the pressure peak occurs, reflected light 
increases, indicating that the layer was dilated during 
the free flight, and that a small compression occurs due 
to the collision. Furthermore, when the layer takes off, 
surface dilation starts to increase (a decrease in reflected 
light). 




0-10 0-12 0-14 0-16 0-18 0-20 

time (s) 



FIG. 9. Time series obtained from pressure and intensity 
sensors for / = 40 Hz; a) F = 1.2, b) F = 2.3. 

We first emphasize that this increase in surface dila- 
tion is the result of the amplification of small differences 
in initial conditions, for the free flight of the grains lo- 
cated at the layer surface |l5|. Such differences arise 
as a consequence of the random character of kinetic- 
energy injection; due to the random packing of grains, 
a layer plate collision naturally induces velocity fluctu- 
ations in the layer. As shown in Fig. || b for F > 2, 
the layer never reaches a compact state, which would 
correspond to a higher value of the reflected light (See 
Fig. ^ a). This result indicates that kinetic energy, in- 
jected into the internal degrees of freedom of the layer, 
has not been completely dissipated within the cycle. 
To estimate the amount of energy not dissipated, we 
consider small changes of / in time for early stages of 
the layer expansion. We can then write I{t)/Iref ~ 



l-2Ads{t)/{d + Ss{0)). RereA6s{t)=6s{t)-ds{0) and 
the reference intensity Iref is taken when the layer starts 
to expand and corresponds to a finite dilation 6s{0) at 
t ~ 0. From the slope of intensity versus time for early 
stages of the layer expansion we obtain a characteristic 
time T. Dimcnsionally t^^ ^ AVto/{d + 6s{Q)), where 
AVto can be associated with the velocity fluctuations at 
the taking off time of the particles located at the free sur- 
face of the layer. Experimentally, r is close to 1/27 s and 
is almost independent of frequency. Thus, our estimate 
shows that AVto/Vc increases linearly with frequency like 
2 X 10-*/, varying from 0.01 for low / (^ 35 Hz) to 0.04 
for intermediate / (~ 200 Hz). Here Vc is the layer-plate 
relative velocity at the collision calculated from the com- 
pletely inelastic ball model. Therefore, our results indi- 
cate that for 2 < F < 2.8, the ratio of residual energy to 
energy injection, (AVto/Vc)^, increases with /, implying 
that energy dissipation decreases with /. At this stage, 
the important feature of typical velocity fluctuations, or 
"temperature", at the free surface of the layer arises: we 
notice that although most of the energy is dissipated, a 
small amount of residual "thermal energy" is enough to 
sustain surface dilation. 

We focus now on the transition from a compact to a 
dilated state suffered by the layer at F ^ 2. Fig. |l^ il- 
lustrates the pressure, P, the collision time, Tc, and the 
reflected intensity versus F for / = 40 Hz. P corresponds 
to the maximum pressure exerted on the plate during the 
collision and Tc here is defined as the width of the pres- 
sure peak at a quarter of its height. We also present in 
Fig. the numerical fit of P with Vc/Tc, it is evident 
that P ^ Vc/Tc, which is a natural consequence of the 
impulsive nature of the periodic forcing. At F ~ 2, the 
DC component of the intensity exhibits a strong decrease 
while its AC component abruptly increases. At the same 
value of F there is a small decrease in the pressure peak, 
and a small increase in Tc, due to an increase in the bulk 
dilation in the layer. This decrease in pressure was al- 
ready observed by P. Umbanhowar |^ , and it is stronger 
for particles with higher restitution coefficients; however, 
no correlation with reflectivity measurements were made 
to investigate the state of the layer. We emphasize that 
the former transition occurs for F < F^,. At the onset of 
surface waves (F = F^,), the pressure presents a strong 
decrease associated with the fact that layer-plate collision 
is spread out in time ||] (Fig. |l^ a, b). 

Using the considerations introduced in the previous 
section, we calculate Sb and Sg as functions of F. This 
is presented in Fig. |ll| for / = 40 Hz. To complete the 
data, we also present the average dilation ((5), and we 
observe that the agreement with Sb is fairly good. Sim- 
ilar to what we found for {6) (fig. P), we observe that 
both Sb and Ss suffer transitions at F w 2. Curiously, 
for F < 2 the bulk dilation is higher than the dilation in 
surface and ds takes negative values, which simply means 
that the layer surface reaches a state more compact than 
the initial one. This is related to the fact that the ini- 
tial state, consistent with our experimental method M, 
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is not the more compact accessible state of the layer. 
Therefore we say that for F < 2 the state of the layer 
is solid-like, where the injected energy during the colli- 
sion is completely dissipated. However, in this regime, 
the available energy is enough to produce rearrangement 
of surface grains. In all the cases tested, the maximum 
compaction was never larger than 2% of the initial den- 
sity. Consequently, we associate the increase in ds and Sb, 
observed at low / and F > 2 to a solid-liquid type transi- 
tion. In the liquid phase, average dilation is large enough 
to allow particles to move with respect to each other. For 
low /, at the critical value F « 2, the injected energy rate 
becomes larger than the dissipation rate, and the energy 
excess sustains the dilation in the granular layer. 



quantities presented in Fig. 10. As opposed to the case 
of low /, at high / and for F > 2 the DC component 
of / increases slightly. This indicates that the layer sur- 
face has reached a state more compact than the initial 
one. It is very important to notice that the decrease in 
pressure associated with the wave instability is clearly 
observed over the entire range of frequencies (Fig. |lj a). 
Complementing Fig. we present in Fig. ^ a and c 
the maximum pressure and the DC component of inten- 
sity versus / for two values of F, both smaller than the 
critical one for waves, with one right below and the other 
right above the fluidization transition. Both the maxi- 
mum pressure and the jump in reflectivity decrease as / 
increases. 




FIG. 10. Maximum pressure P, collision time Tc, DC and 
AC components of intensity versus F for / = 40 Hz. In (a) 
we present the fitting of the maximum pressure P (•) with 
Vc/Tc (A), whose magnitude is presented in the right axis in 
units of m/s^. 
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FIG. 11. Bulk dilation, Sb (•), mean dilation ,{S) (A), and 
surface dilation, Ss (*), versus F for / = 40 Hz. 

In Fig. Owe include frequency dependence of the same 




FIG. 12. Maximum pressure, P, collision time, Tc, and DC 
and AC components of intensity versus F for various f;f = 
40 (•), 69 (A), 83 (x). 111 (o), 154 (+), 200 (o), 250 (*) and 
350 (□) Hz. 

We also present the bulk dilation Sb (See Fig. |l^ b) 
calculated directly from the pressure through equation 
d). We find that St scales as with 6 = -1.42 ± 0.07 
and b = —1.54 ± 0.03 for F = 1.5 and 2.3 respectively. 
These values are in very good agreement with those ob- 
tained for (i5) (See Fig. ^), and they provide a consis- 
tency test for both kinds of measurements. These results 
indicate that the relevant quantity is not only the ratio 
of the injected energy per particle to the potential en- 
ergy required to rise a particle by a fraction of its diam- 
eter, in which case Sb would vary as l/f^, but also the 
dissipated energy. We notice that our results contrast 
with those obtained by Luding et al in numerical simula- 
tions of a column of particles in the completely fluidized 
regime, where the average dilation scales as {5)j;^ (Af)'^, 
which for F constant becomes (S) ~ {T / How- 
ever, similar simulations done by the same group for a 
two-dimensional layer |l^ indicate that the layer expan- 
sion scales as hem — hcmo (Af)^^^ ^ (F//)^/^. Even 
though these results were also obtained in a completely 
fluidized regime (typically F > 10) we observe quite good 
agreement for the scaling on f for the expansion of the 
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layer. 

On the other hand, information about the surface di- 
lation, Ss, versus / is obtained using the intensity data 
as I/Iq « (f/id + Ss f (See Fig. |l| d). For / > 225 
Hz, Ss takes negative values, which indicates that the 
surface layer is more compact than the initial one [Q. 
However, in this regime, as shown in the previous sec- 
tion, the bulk dilation increases with F but remains very 
small, d/d ~ IQ-^. 
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FIG. 13. Maximum pressure P, bulk dilation St, DC com- 
ponent of intensity and surface dilation versus / for two 
constant values of F below the onset of surface waves; F 
= 1.5 (•), 2.3 (A). Bulk dilation scales as 5b ~ Z*" with 
b = -1.42 ± 0.07 (•) and b = -1.54 ± 0.03 (A) respectively. 
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FIG. 14. Maximum pressure P (a) and DC reflectiv- 
ity-component (b) as a function of F for / = 60 Hz. The 
vertical line indicates the flat layer with kinks transition at 
which Vc = 0. At F = 3.6, as predicted by the ineslastic col- 
liding ball model for the center of mass of the layer, a period 
doubling instability is observed. 



Finally, let us mention another interesting transition 
which links to the flat with kinks instability reported in 
previous works |Q . If we increase F further we find that 
a period doubling is achieved at F w 3.6, as reported 
before ]l7|,pt. Next, for F « 4.6, an inverse transition of 
the liquid-solid type is detected. Fig. ^ shows a typ- 
ical measurement of both reflected intensity and maxi- 
mum pressure for a wider range of F. We detect that 
at F « 4.6, P strongly decreases and the mean value of 
intensity strongly increases. Using the results discussed 
previously, both changes reflect a strong decrease of grain 
mobility. We notice that the critical value F w 4.6 corre- 
sponds to Vc = in the completely inelastic ball model, 
so in fact no energy is injected to the internal degrees of 
freedom. We also conclude that no surface waves can be 
sustained in this regime, except at the kink itself. Indeed, 
the shear induced at the kink by the flat parts oscillating 
out of phase is large and enough to induce dilation. At 
low frequencies, this dilation is enough to allow hydro- 
dynamics waves, as those shown in figure As F is 
increased further, the pressure increases and the reflec- 
tivity decreases. Thus, energy injection again becomes 
enough to sustain surface waves in the layer. This fact 
is consistent with the existence of //4 waves reported 
previously |^]. Notice that for F = 4.6 the maximum ve- 
locity of the plate is Aw « 12 cm/s; this confirms that 
the relevant scale of velocity fluctuations is given by Vc 
and not Aw. 

The experimental results presented above can be sum- 
marized as follows. Depending on the excitation fre- 
quency we observe different kind of states and waves. At 
low frequency, bulk and surface dilation present strong 
increases which are associated with a fluidization tran- 
sition. Surface waves observed in this regime involve 
large relative motion between particles and are there- 
fore considered as the hydrodynamic modes of the layer. 
At intermediate /, although the injected energy is small, 
we still observe a decrease in reflectivity as a function 
of F. In this regime, as shown in Fig. |l^ b, St is too 
small (db/d < 0.1) to allow motion between the parti- 
cles [Q. Therefore, at the critical F, the decrease of re- 
flectivity is the signature of particles fluctuating around 
their positions at the free surface of the layer. We asso- 
ciate this decrease with a heating up of the solid phase. 
For higher frequencies, the layer undergoes a compaction 
transition which is detected by the increase in surface 
density. Below and above this transition the local bulk 
dilation Sb/d is even smaller and is of order 0.005, im- 
plying that the mobility in both the bulk and the free 
surface is completely suppressed. Thus, very low ampli- 
tude surface waves detected in the compaction regime, by 
the strong decrease in the maximum pressure, must cor- 
respond to excitations in which the layer is slightly mod- 
ulated in time and space. We will see in the following 
section that these waves are bending waves, associated 
with the ability of the compact layer to deform. Finally, 
fig. |l^ presents the phase diagram for the granular layer: 
phase boundaries separating the various states and sur- 
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face waves have been obtained from the data in Fig. 
The layer state and surface wave transitions occur for 
approximately constant values of F independent of /. 
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FIG. 15. Phase diagram showing layer state and surface 
wave transitions. The vertical dashed line at / — 225 Hz de- 
fines the frequency above which the layer strictly undergoes 
a transition to a state more compact at F ~ 2. 



V. LOW AMPLITUDE WAVES: BENDING 
WAVES 

To check the existence of these waves we have per- 
formed a set of experiments in a two-dimensional gran- 
ular layer. The advantage of using a two-dimensional 
system is that this allows us to obtain side views of the 
waves. This is of particular importance since bending 
waves are difficult to visualize as they have amplitudes 
of the order of fifty percent of a particle diameter. 

In order to observe such small amplitudes we have used 
photoelastic cylinders of 6 mm in diameter and 6.35 mm 
in length. Typical excitation frequencies arc about 40 Hz, 
which for this system are in the high frequency regime. 
We recall that this regime occurs for frequencies much 
larger than the frequency crossover, / ~ y/gjd. When 
varying the deepness of the layer (number of particles = 
N), this scaling becomes / ^ \fgJWd [|l9|. In the case 
of bronze particles of d « 0.12 mm, bending waves are 
detected for / > 225 Hz. This tells us that these waves 
should be observed at / « 40 Hz, for d = 6 mm and 
a layer about ten particle diameters thick. With these 
parameters, for F ^ 3.5, the amplitude of waves should 
be of the order of 1 mm. 

Typical snapshots of two stages of bending waves at 
F = 3.5, / = 40 Hz and iV = 10 are presented in figures 
[T6| a and p^b. We observe that the layer slightly bends 
with respect to the horizontal. Similar to what occurs 
for low frequency waves, this modulation alternates in 
time at half of the frequency forcing (See below). How- 
ever, in this case, the wavelength of the modulation is 
about a layer thickness and is also nearly independent 
of /. Some particles are marked with black spots which 
allows us to distinguish them and follow their trajecto- 



ries (See Fig. [l^a, b). As expected, the mobility of the 
particles is very low. Only at the surface layer will some 
particles move with respect to each other over distances 
of the order of the driving amplitude. In the bulk this 
motion is completely suppressed. 




FIG. 16. Typical snapshots of bending waves for / = 40 Hz 
and F = 3.6. In (a) and (b) the white horizontal fines show 
that the layer bends itself about half a diameter. Both states 
alternate in time with the frequency of oscillation. In (c) we 
present the difi'erence of two images during a collision. The 
bright zones correspond to regions under high stress. In this 
case the compression front travels to the left. 

Additional support for these low amplitude waves is 
provided in Fig. |l6| c in which a compaction front that 
moves laterally is observed as a bright zone. This image 
results from the difference of two consecutive snapshots 
(period of acquisition is about 0.8 ms). At a compression 
zone, the vertical stress in the layer is high so the light is 
transmitted. Thus, the difference of images mainly shows 
the zones under stress. In the case presented the com- 
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pression front moves from right to left and the coUision 
occurs very near the right boundary of the image. In 
the next cycle the collision occurs near the left boundary 
and the compression front travels to the right. With a 
wider view of the layer we observe that in fact two com- 
pression fronts are created from each collision point; one 
front travels in each direction. This kind of visualization 
allows us to safely say that these parametric waves are 
also subharmonic, i.e. their frequency is //2. Finally, 
from the images we estimate the velocity of the compres- 
sion front of the order of a few meters per second. This 
value is very high compared to the estimated velocity at 
the collision Vc ~ 0.25 m/s; this indicates the existence 
of contact arcs (see the contact lines in Fig. |l^) so the 
layer is almost not dilated. 



VI. CONCLUSIONS 



In conclusion, depending on the excitation frequency, 
we observe different kinds of states and waves. Thus, 
our experimental results reveal the existence of a solid- 
liquid transition that precedes subharmonic wave insta- 
bility. Hydrodynamic surface waves can be then con- 
sidered as the natural excitations existing in a fluidized 
granular layer. In contrast, very low amplitude surface 
waves detected in the compaction regime correspond to 
excitations in which the layer slightly bends alternatively 
in time and space. We have seen in the previous section 
that these waves are associated with the compact charac- 
ter of the layer. The layer states observed here and sur- 
face wave transitions occur for approximately constant 
values of F independent on /. 

Additional experimental evidence of the fluidization 
transition experienced by the layer at F ~ 2 can be 
found in reference (See for instance, Fig. 6 in ref. 
pO|). Although they used rather large particles at small 
forcing frequencies, we find that their experiments are 
actually in the low frequency regime after estimating the 
frequency crossover. Independent experimental evidence 
for the compaction transition is also found in previous 
works. For instance, a strong increase in granular den- 
sity was found, close to F ~ 1.9, in several experiments 
on large columns of grains submitted to "taps" of a sin- 
gle cycle of vibration |21|. In this case, applying the 
scaling to the frequency crossover, we found consistently 
that these experiments correspond to the high frequency 
limit. Therefore, we conclude that both fluidization as 
well as compaction transitions are well established ex- 
perimentally. However, it is still unclear what mecha- 
nisms dominate these transitions and why they arise at 
a constant value of F. We can only safely say that at 
low energy injection rates, or equivalently at small plate 
amplitudes with respect to the thickness of the layer, the 
transition will be of the compaction type. In the opposite 
case of high energy injections rates, this transition will 
be of the fluidization type. 



Also, we have clearly shown that the relevant scale of 
velocity fluctuations is given by Vc and not Aw. Then, 
via the intensity measurements we deduced that the dis- 
sipation decreases as / is increased; this is equivalent to 
state that the dissipative effects increase with velocity 
fluctuations. A possible cause of this is the dependence 
on the velocity of the restitution coefficient: it has been 
recently shown that 1 — e ~ ti", where v is the relative 
normal velocity at the colision and a a positive number 
|p2| . Nevertheless, as the dilation is reduced we expect 
that the friction between grains will become an impor- 
tant dissipative mechanism. 

Finally, the importance of the bulk and surface dila- 
tion measurements is that they provide a complementary 
way, with respect to the granular temperature, to explore 
the excitation of the internal degrees of freedom of a vi- 
brated layer. For the dilation in the bulk we have found 
that, for a constant value of F, it is a decreasing function 
of / of the form 1//'', with b « 3/2. This agrees with 
previous simulations of a two dimensional layer p^ ] . It is 
clear that the deviation from the expected value 6 = 2 is 
due to dissipative effects, but the exact numerical value 
seems to depend on F. 
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